Supercell band calculations and correlation for high- Tc copper oxide superconductors. 
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First principle band calculations based on local versions of density functional theory (DFT), 
together with results from nearly free-electron models, can describe many typical but unusual prop- 
erties of the high- Tc copper oxides. The methods and a few of the most important results are 
reviewed. Some additional calculations are presented and the problems with the commonly used 
approximate versions of DFT for oxides are discussed with a few ideas for corrections. It is con- 
cluded that rather modest corrections to the approximate DFT, without particular assumptions 
about strong correlation, can push the ground state towards anti-ferro magnetic (AFM) order. Spin 
fluctuations interacting with phonons are crucial for the mechanism of superconductivity in this 
scenario. 

PACS numbers: 74.25. Jb,74.20.-z,74.20.Mn,74.72,-h 



I. INTRODUCTION. 

It is often assumed that strong correlation is irnportant 
for an understanding of the high-Tc problem [l|, |3, S|- 
The failure of the local approximations to the exchange- 
correlation functional, such as the local (spin) density 
approximation, LDA or LSDA [1,0| to produce the 
anti-ferro magnetic (AFM) insulating state of many un- 
doped transition metal oxides is generally quoted as the 
essential reason for discarding DFT calculations for high- 
Tc copper oxides [3j. DFT is essentially exact, but the 
approximations to make it practical for use in real appli- 
cations seem inappropriate for oxides. However, despite 
this problem there are several LDA results that fit to the 
observed high-Tc properties, provided that doping and 
supercells are considered in order to account for imper- 
fect lattice conditions such as stripes, phonons or spin 
waves [1, ■ Here is presented a short review of those 
DFT results together with a discussion of the correlation 
problem and an attempt to include further corrections to 
LSDA. Estimations of couplings A caused by spin-phonon 
coupling (SPC) and pure spin fluctuations are also made. 



II. NEARLY FREE-ELECTRON MODEL. 

The spin polarized potential from an AFM order 
(like that of undoped La2Cu04) can be generated from 
VAFMix) = Voexp{—iQ ■ x), where Q = tt/oo, is at the 
Brillouin Zone (BZ) boundary, so that there is one Cu 
site with spin up and one with spin down within a dis- 
tance ao (ao is the lattice constant). The potential of 
the other spin has a phase shift of tt. The band dis- 
persion from a 1-dimensional (ID) nearly free-electron 
model (NFE) is obtained from a 2x2 eigenvalue prob- 
lem, and it has a gap of 2Vo at the zone boundary. 
Suppose now that an additional potential modulation 
Vq{x) = Vqexp{iq - x) exists, and that | q\<^\ Q |, so that 
its periodicity or "wavelength" is much larger than 2ao. 
The product of these two modulations can be written 
V{x) = Vqexp{—i{Q ~ q) ■ x), where the two amplitudes 
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FIG. 1: A picture of how stripe like configurations are made 
up from the product of two plane wave potentials. The arrows 
in the upper row show the AFM moments on Cu sites given 
by the envelope function Re{exp{—iQ ■ x)) (broken line), and 
the unit cell contains 2 Cu sites along x. The second row 
shows the modulation Re{exp{iq ■ x) (here q = Q/4) and the 
last row is the product Re{exp{—i{Q — q) ■ x) with the new 
spin configurations on the Cu. The striped unit cell contains 
8 Cu sites. 



are combined into one coefficient, Vq. Fig. 1 shows the 
real space configurations along the 100-direction. This 
potential describes AFM modulations in stripes, with 
magnetic nodes (Cu with zero moment) separated by 
aoQ/q. The values of Vq for phonons or spin waves of 
different lengths are not known from the NFE model, 
but some values will be fed in from the ab-initio Linear 
Muffin-Tin Orbital (LMTO, [ll|) calculations, as wiU be 
described below. The LMTO calculations are slow for 
large supercells. The cells need to be large to cover the 
periodicity of realistic phonon and/or spin- waves, and so 
far it has been possible to extend the cells in one direction 
only, usually the CuO bond direction. In contrast, the 
NFE model is very simple and it can easily be extended 
to two dimensions (2D). The band dispersion for k-points 
{kx, ky) is now obtained from the eigenvalues of a 3x3 ma- 
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FIG. 2: Examples of the DOS from the 2D NFE model for 3 
different strengths of Vq (in mRy) for (0.075, 0.075). 



trix Not only stripe order, but also "checkerboard" 
configuration can be modeled. 

Examples of the density-of-states (DOS) for 3 different 
Vq are shown in figure 2. The potential modulations are 
parallel to the CuO bond directions along x and y, so 
the gaps appear on the {k^, 0) and (0, ky) lines, as shown 
in figure 3. However, almost no effect from Vq appears 
on the band dispersion in the diagonal direction (along 
{k, k)). Therefore, the Fermi surface (FS) is not affected 
along the diagonal, but fluctuations, described through 
amplitude variations of Vq for different positions and dif- 
ferent time, will make the FS smeared in the two bond 
directions. An example of this is shown in figure 4. 



III. AB-INITIO ID-LMTO. 



FIG. 3: An example of the 2D NFE bands along (fe, 0) and 
{k,k) for qx = qy = 0.075 and Vg = 20 mRy (full line). The 
FE band without potential perturbations {Vq —0) is shown 
by the broken line. The thin horizontal line is at the DOS 
minimum, which is situated at 0.15 holes/cell. 



By ID (1-dimensional) LMTO we mean ab-initio 
LMTO calculations for long (and narrow) supercells most 
often oriented along the CuO bond direction [l^ 
These calculations are based on the local version of DFT, 
the local (spin) density approximation, LDA or LSDA 
0, [^. Cells with phonon distortions and/or spin waves 
within lengths of 4,8 and 12 lattice constants are typi- 
cally considered in these calculations. The wavelengths 
of spin waves are twice as long as those of phonons. Hole 
doping, ft,, in La(2-h)BahCu04 (LBCO) is generally mod- 
eled by the virtual crystal approximation (VCA) where 
the nuclear and electronic La-charges (57.0) are reduced 
to {57-h/2) to account for a perfectly delocalized dop- 
ing {h in holes per Cu). But some calculations with 
real La/Ba substitutions show that improved supercon- 
ducting properties can be expected from periodic dop- 
ing distributions p^ . The maximal phonon distortion, 
Ui, for different sites, i, and AFM magnetic moments on 
Cu, m, depend on temperature, T, force constants and 
spin stiffness. Appropriate values of Ui{T) and m(T) at 
T w lOOK are deduced from experiments and calculations 



FIG. 4: An example of the Fermi surface in the 2D NFE 
model with fluctuations of Vq. The FS is well defined in the 
diagonal direction, while the fiuctuations makes it diffuse in 
the X and y-directions. 



[H, [BEES mi. 

A striking result of these calculations (also made for 
HgBa2Cu04) is that a gap (or a pseudogap) will open 
up in the DOS, very similar to what is found in the NFE 
models. The gap will appear near to Ep for the undoped 
material if the wavelengths (A) are long, while for short 
phonons or spin waves the gap moves to lower energy 0] . 
In summary, A = 1/h, where A is in units of oq and h is 
the number of holes per Cu. 

A second important result is that the spin waves will 
be stronger, and the pseudogaps will be deeper, if the 
waves coexist with phonon distortions of the correct wave 
length and phase [Tsj . Interactions between phonons 
and spin waves have also been suggested in order to ex- 
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plain neutron scattering [13, H^l ■ Phonons like the " half- 
breathing" 0-mode and modes with z-displacements of 
La and apical oxygens are most effective in this pro- 
cess of SPC. This mechanism offers direct explanations 
of phonon softening, q-dependent spin excitations and 
various isotope effects [lO|, [ij] . 

The rather few LMTO results permit to establish only 
a few values of Vq for different q,i of the phonons and 
spin waves. A general trend of larger Vq for long wave 
lengths seems clear, both for phonons and for spin waves. 
A procedure based on the partial character of the states 
above and below the gap in the NFE model confirms 
this trend for spin waves, and it permits to fill in some 
of the q-dependent points not calculated by LMTO. In 
all this gives some confidence in the q, i-variations of Vq. 
However, one can expect more vivid variations of the 
spin part near g — > 0, because LDA underestimates the 
transition towards AFM. 

When these Vq-values are fed into the 2D NFE model 
it gives an approximately linear variation of q as function 
of h, up to a saturation near q ~ 0.125 for h larger than 
^ 0.13 [13 . Hence, the linear dependence for low h is in 
qualitative agreement with the result from ID-LMTO, 
but the pace is different. This is probably because the 
gap opens both along x and y in the 2D NFE model, 
which makes the progression of the gap a bit slower than 
in ID. The saturation is because the Vq^s decrease for 
increasing q and it becomes impossible to open a gap at 
the low energy where Ep should be for large h. How- 
ever, if qx and qy are assumed to be different in the 2D 
NFE model it is possible to follow one gap towards lower 
energy. A second weaker gap moves to higher energy. 




It was suggested that non-adiabatic electron-phonon 
coupling could be enhanced in the cuprates, because of 
the low Fermi velocity in the z-direction 24]. This ve- 
locity is comparable to the vibrational velocity and the 
electronic screening appeared to be insufficient. How- 
ever, the screening can be made by other electrons mov- 
ing within the planes. Moreover, phonon frequencies cal- 
culated within LDA without assumptions of incomplete 
screenin g, agree satisfactory with experimental frequen- 
cies [l8lTl9l.l20|. Instead, excitations of virtual phonons 
coupled to spin waves can be important for the mecha- 
nism of superconductivity, but SPC makes the separa- 
tion between pure electron-phonon coupling, Aep, and A 
caused by spin fluctuations, A^/, less clear. 

The important observation is that atomic phonon dis- 
tortions will trigger enhancement of spin waves. If so, 
the possibility for larger Agp is open, because instead for 
the common approach to ignore spin effects in electron- 
phonon coupling, there are larger matrix elements when 
the spin polarized part of the potential is involved. This 
is readily imagined for a system which is nonmagnetic 
when phonon distortions are absent, but magnetic when 
the distortions are present. An estimation of pure Aep 
and the coupling parameter for SPC, Xspc, in LBCO 
has been presented earlier With a total N{Ep) « 
0.9{eV ■ cell ■ spin)~^, and distortion amplitudes and po- 
tential shifts as in ref. [l2| this leads to X's of the order 
0.36 and 0.6 for pure phonons and SPC, respectively [2^ . 
These values appear sufficiently large for a large Tc in 
simple BCS-type formulations, but the precise relation 
for Tc depends also on other parameters [1^ . 



IV. SUPERCONDUCTIVITY 

The emerging picture from these band results is that 
moderately strong spin fluctuations, which exist for many 
combinations of q and uj, will be enhanced through SPC 
to some particular phonon distortions. The selection de- 
pends on the atomic character of the phonon mode as 
well as on doping and wave length (i.e. on q). The lat- 
ter is because the pseudogap appears at Ep so that a 
maximum amount of kinetic energy can be gained from 
the SPC mode. This is for low T when the states be- 
low the gap is of one spin and well separated from the 
(unoccupied) state of the other spin above the gap. But 
for increasing T there will be mixed occupations of the 
states around Ep through the Fermi-Dirac function (and 
through thermal disorder), which will decrease the spin 
density. This will decrease the spin polarization of the 
potential, which in turn will decrease the spin density 
even more, and at some temperature T* the support of 
the pseudo gap from the spin wave will collapse [l^] ■ The 
fact that Vq is largest at low doping favors larger T* when 
h ^ 0. A high DOS at Ep is important for a high super- 
conducting Tc, and therefore is the pseudogap in com- 
petition with superconductivity in under-doped cuprates 



TABLE I: Total energy {AE, mRy per Cu), matrix element 
for spin fluctuations {AV in mRy), As/ and local Stoner 
enhancement, S, calculated for lattices of LBCO with two 
phonon distortions and without distortions. 



phonon 


AE 


AY X,f 


S 


no phonon 


12.7 


9 0.03 


1.6 


plane-O 


8.2 


12.5 0.12 


2.6 


La 


7.4 


12.0 0.11 


2.4 



The third mechanism is that spin fluctuations work 
without coupling to phonon excitations. Still, phonon 
distortions might be present, and will in that case be 
important for enhancements of rapid (high energy) spin 
fluctuations. For instance, it was found that the local ex- 
change enhancement (for AFM) on Cu sites are largest 
when the surrounding atoms (La or oxygens) have been 
pushed away from the Cu, and this leads to SPC be- 
tween phonons and spin waves with equal q and to [l^. 
But strong high-w spin excitations should also be possible 
on rows of Cu with enhanced exchange, perpendicular to 
the phonon. Ab-initio calculations are not easy for such 
configurations because of the large size of the required 
unit cells, but one can use some of the existing results for 
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large w. In contrast, SPC is essentially limited to below 
the highest phonon frequency, i.e. to 50-70 meV. The rel- 
atively small absolute values of Xsf are partly caused by 
the use of LDA. Improved DF schemes which could pre- 
dict an AFM insulating state for the undoped systems, 
would lead to larger exchange enhancements, larger A^/ 
and lower Wg/ in doped systems. 



CORRELATION AND CORRECTIONS TO 
LDA. 



to (meV) 

FIG. 5: The relative strength of Xspc (the first 4 points) 
and As/ (the 3 points at the highest energies). The latter 
will move towards lower energy if the exchange enhancement 
increases, while the points for SPC are fixed near the phonon 
energies, see the text. 



first estimations of Xsf for rapid spin fluctuation (so rapid 
that the phonon distortions appear to be static compared 
to the spin fluctuation). The enhancements depend on 
the different type of distortions as was discussed above. 
In this case we calculate Xgj — N{dV/dm)'^ /{(PE/dm?), 
where V is the potential, E is the total energy of the 
spin wave and m is the magnetic moment (per Cu) [26| . 
The difference in free energy between non-polarized and 
an induced (by magnetic field) AFM wave can be writ- 
ten E,n ~ Eq + Kui^, so that d^E/dm^ — 2k, which 
permits to calculate Xgf — N ■ {AV)'^ /k for "harmonic" 
spin fluctuations. With parameters from the band results 
this makes Xsf equal to 0.03 when no lattice distortions 
are present, and 0.12 and 0.11 when the lattice contain 
distortions, see Table 1. 

The frequency has not been calculated explicitly, but 
a shortcut via exchange enhancements can be used for a 
simple estimation, as for FM fluctuations [l^ll^]. Thus, 
ujsf « l/(4iVS'), where S = AV/AH is the local Stoner 
enhancement (on Cu) and AH is the external field ( 5 
mRy in these calculations). This makes uisf 100-200 meV 
depending on the type of lattice distortion. Fig. 5 shows 
the combined results for direct Xspc and indirect As/ as 
function of frequency. The general shape is similar as the 
electron-boson coupling function that has been extracted 
from recent optical spectra on Hg- and Bi-based copper 
oxides 01 . 

The amplitudes of Xspc and Xsf increase further if 
both of the latter distortions (plane oxygens and La) are 
considered in the calculation. The matrix element in- 
creases to about 14.5 mRy. Despite this relative increase 
it is seen that pure spin fluctuations give not as large X's 
as for SPC. This is especially clear when all lattice modes 
are considered. However, these pure spin fluctuations can 
be important for superconductivity, since they appear at 



The question is if the real electron-electron correlation 
in the oxides is much larger than what is included in 
LDA. The latter includes exchange and correlation (XC) 
for an electron gas of varying density with the electron 
gas radius = 0.62p~^/^ as parameter. The density p 
contains one electron within rg. The largest values of 
(« 2a. M.) are found in the low density valence-electron 
region between the atoms, but they are still smaller than 
any typical atomic radius, Ra in transition metals, their 
oxides and these cuprates. This gives a major argument 
against using a strong on-site correlation parameter to 
an atom, because the effect of an additional electron will 
readily be screened out, and is not noticed beyond rs- In 
low-density materials however, such as the alkali metals, 
rs ~ Ra and on-site correlation could be advocated. But 
LDA seems to describe alkali metals well, which provides 
an additional argument in favor of local approximations 
of DPT. 

The fact is that LSDA does not produce an AFM in- 
sulating ground state of several oxides, so the problem is 
real and it may be related to the spin-polarized part of 
the potential. Small errors with important consequences 
of this type are known for LSDA. A well-known example 
is for Fe, where LSDA predicts nonmagnetic (NM) face- 
centered cubic (fee) lattice as the ground state, whereas 
the generalized gradient approximation (GGA, ^]) cor- 
rectly finds the ferro rnagnetic (FM) body-centered cu- 
bic (bee) ground state [2^ . The difference between LDA 
and GGA is not large concerning the bands, magnetic 
moments and other properties. The band properties of 
bcc Fe are quite good in calculations using both these 
DFT functionals. But the order of the total energies for 
FM-bcc and NM-fcc is reversed in the two cases. 

An additional indication that the AFM insulating state 
of the cuprates is not far away in LSDA is the fact that 
the AFM moment depends on the number of k-points. A 
sufficient number of k-points is required for convergence, 
but early works noticed that the AFM ground state be- 
came stable if a coarse k-point mesh was used [13, HH . A 
similar conclusion for the appearance of weak ferromag- 
netism in highly doped LBCO can be understood from 
DOS effects and the Stoner criterion [11]. Thus, as for Fe, 
a small detail can have large consequences. It might be 
that the total energy difference between two very differ- 
ent ground states (fee nonmagnetic and bcc-magnetic for 
Fe, and NM metallic and AFM insulating for cuprates) 
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are very close, so that small errors in the computational 
details lead to the 'wrong' ground state. 

Another indication of a nearby AFM state is that 
LMTO calculations with off-center linearization energies 
can stabilize AFM in HgBa2Cu04 in calculations with 
sufficient number of k-points. The normal procedure in 
LMTO is to chose linearization energies at the center of 
the occupied sub bands, but if they are chosen at the bot- 
tom of the bands, it leads to a slight localization of the 
states. An AFM order needs less hybridization energy 
and AFM can be stabilized [s^. This is rather ad- hoc, 
since LMTO should not be made in such a way. A theo- 
retical justification is need for a better DFT scheme. 

A first workable non-local density functional is the 
GGA (generalized gradient correction Q), where the 
gradient, or the first derivative of the charge density, 
dp{r)/dr, brings out corrections to the LSDA. The gradi- 
ents are mainly radial in atoms and solids because of the 
steep drop in charge density for increasing r close to an 
atom. Only valence electrons have large amplitudes in 
the interstitial region between the atoms, where the gra- 
dients are moderate. However, higher gradients such as 
the second derivative of the density can be large in the 
interstitial, a feature not captured by the GGA. Nev- 
ertheless, by using GGA in the LMTO calculations for 
LB CO there is some improvement, since the required crit- 
ical magnetic field for having a gap is reduced by about 
25 percent compared to LSDA. 

Densities with quadratic dependence as function of r 
around a fixed point at r = are used in a 2-particle 
model for calculations of correlation including corrections 
due to the second derivatives. A non- interacting density 
with quadratic gradients have the form 

Pi^)=Po + ^^V^ (1) 

The Schrodinger equation is solved for a density of elec- 
trons surrounding a fixed electron at r = 0, where it is 
taken into account that the effective mass is 1/2 [33.[35|. 
The interacting potential can be written 

Vir) = eVr + £{£ + l)/r^ + Ai,,(r) + K.t W (2) 

The exchange-correlation within the surrounding elec- 
tron cloud is taken into account through LDA (fj-xc), and 
the external potential Vext can include relative kinetic en- 
ergy variations. The two last terms in V{r) are not very 
important for the results for gradient variations. The 
strongest term is the unscreened Coulomb repulsion that 
diverges at r = 0. The second term, the centrifugal term 
is included only for exchange (with i = 1, higher £ will 
have higher energy), when the Pauli principle requires 
that the density for equal spin must be zero at r = 

The solution ^'(£^,r) of the Schrodinger equation for 
this potential is used to calculate the XC energy as the 
difference in Coulomb energy between interacting and 
non-inteacting densities; fi = /(^'(£',r)^ — p{r))/rd^r. 



The energy E is determined from a required boundary 
condition of '^{E,rs). The variation of the correlation 
potential on the second gradient is quite easy to under- 
stand. If the gradient is positive, so that there is a ten- 
dency for a 'hole' at r—Q, it will be less effective for the 
Coulomb repulsion to make the correlation hole deeper at 
this point. Therefore, fic will decrease for positive density 
gradients near the interstitial region. The gradient can 
change closer to the high density regions near the nuclei 
and make the correlation larger than for a constant den- 
sity, but correlation becomes negligible in comparison to 
exchange when the density is high. The C-correction is 
parametrized through the electron gas parameter and 
Q = p{r,)/p{0), so that C = 1-5^^/4+ y^*Q/2 for 
Q > I and as the inverse of this parametrization if Q < 1 . 
The expression is only applied for Q between -0.25 and 
0.8, and C — 1 is used as a scaling factor of the correla- 
tion part from LDA. This correction makes the potential 
a little more repulsive in the interstitial region. The Cu-d 
states becomes more localized, which can promote AFM 
order, see Table IL 

One fundamental theorem of DFT is that of "v- 
representability" , that essentially says that there is a 
one-to-one correspondence between the density p{r) and 
potential V{r) Q. The only possibility is that two equiv- 
alent charge densities can give an uninteresting shift in 
the potential, which can be absorbed by the shift of the 
energy spectrum. Imagine now that an external potential 
is applied to an electron gas of varying density p(r). The 
free electron wave functions will oscillate more if more 
kinetic energy is given to the electron gas, or less if ki- 
netic energy is subtracted, even if p remains constant. 
The exchange part of the potential (X) will be modified 
because of the ability of the wave function to readjust 
itself around a second electron of the same spin (to cre- 
ate an exchange "hole"), as can be imagined from Slater 
exchange Q or the two-particle model [s^. A similar 
change occurs for correlation (C), which is caused by the 
immediate Coulomb repulsion between all electron pairs. 
If p is constant in space there is only an uninteresting po- 
tential shift, but for varying charge densities there can be 
more subtle effects. The LDA is derived from the prin- 
ciple of density variations at the Fermi energy, and the 
relation between Ep (which is a measure of the kinetic 
energy) and density is given by; 

£;«(r) = (37r2p(r))2/3 (3) 

(2) 

This relation can be compared to Ep (r) = Ep — V{r) 
within real atoms, where Ep is from the band calculation 
with potential V{r). Except for r very close to the nuclei 

there is not too large a difference between Ep\r) and 

Ep (r) . In some regions the two values can be quite dif- 
ferent. If so, it implies that the density is not in equilib- 
rium with the kinetic energy as was used in the derivation 
of the LDA potential. Oxygen sites in oxides are often 
negatively charged, which suggests that the kinetic en- 
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FIG. 6: The calculated values of the scaling function X as 
function of Vs and /e (see text). 

ergy should be smaller than what comes out from eq. 1. 
These effects are moderated by the fact that the charge 
transfer is a result of an attractive potential. It turns out 
that positive and negative deviations from eq. 1 exist in 
shells within all atoms. 

Thus, band calculations with correction for the kinetic 
energy will be attempted where the local exchange po- 
tential is written 

f^xifcTs) = X{fe,rs)nKs{rs) (4) 

where X{fe, r^) is a scaling function of the normal Kohn- 
Sham potential fiRS, 

A derivation of the scaling function X{fe,rs) can be 
done directly from the Slater functions by use of energy 
shifts in the arguments of the plane waves Q- How- 
ever, negative shifts leading to localized waves have to be 
avoided, and the resulting X-function seems too sensitive 
to small variations of f^. As an alternative we determine 
the scaling function from a two-particle model as in ref. 
[H, UHli with a renormalization to make X(l,rs) = 1. 
The result is displayed in fig. 6 for the appropriate range 
of rg and /g. The real value of /g for r < Q.QbRnws 
increases and can be larger than 50 near the nuclei. This 
is extreme and fe is here cut off at 3. 

TABLE II: The required applied magnetic fields (in mRy on 
Cu) in order to obtain a zero gap in LMTO calculations which 
include the different corrections to the potential, see the text. 



LSDA GGA GGK+C GGA+X GGA+C+X 
8^8 6J 5l) 4l) 



In preliminary calculations for undoped LBCO we do 
a rescaling of the exchange due to kinetic energy and of 
correlation due to second gradients. The comparison is 
made with standard LSDA for an AFM unit cell where a 
staggered magnetic field is applied to the Cu sites, and it 
is found that a field of ± 8.8 mRy is the limit for having 
a zero gap {Eg 0.25mi?y). The magnetic moment is 
then ±0.18/x_B per Cu. The amplitude of the required 
magnetic field to obtain a zero gap is reduced when the 
calculations include the correction factors, see Table II. 
The absolute values can depend on details of the band 
calculations |33|], but the trend towards a stability of an 
AFM state is clear. Corrections of this type will be in- 
teresting for seeing enhanced spin fluctuations for long 
wave length spin waves in doped cuprates. Further en- 
hancements of Xgf can also be expected |l2|. Potential 
corrections must ultimately be tested for other types of 
materials in order to verify that some properties will not 
deteriorate. Here for the cuprates the corrections permit 
to proceed in the modeling of more properties from the 
band results. 



VI. CONCLUSION 

Results of band calculations for supercells with frozen 
phonons and spin waves suggest that a Fermi-liquid state 
can cause pseudo gaps and dynamic stripes. Together 
with a NFE model it is possible to describe the dop- 
ing dependence of many normal state properties of the 
cuprates. The coupling between phonon distortions and 
spin fluctuations seems to be crucial for the mechanism 
of superconductivity, so that the spin-polarized part of A 
is most enhanced by simultaneous excitations of phonons 
and spin waves. Two different mechanisms for supercon- 
ductivity mediated by spin fluctuations are possible. The 
largest coupling parameter is when a phonon is excited 
together with the spin fluctuation. Lower couplings. As/ 
at larger energies, are independent of the phonon excita- 
tion, but these spin fluctuations can nevertheless profit 
from possible phonon distortions of the lattice. Still, ab- 
solute numbers are too small when using LDA, as is also 
concluded from the absence of AFM stability in LDA cal- 
culations for undoped systems. However, it is argued that 
corrections due to higher order density gradients and ki- 
netic energy are able to bring the band calculations closer 
to AFM. This is shown in calculations for undoped LBCO 
by yet very approximate corrections due to non locality 
and kinetic energy. Refinements of such corrections will 
be of interest for application to supercells, including dop- 
ing, phonon distortions and spin waves, since it can be 
expected that realistic A's will be obtained. 
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